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Abstract. Using simple particle models of limit order markets, we argue that 
the mid-term over-diffusive price behaviour is inherent to the very nature of these 
markets. Several rules for rate changes are considered. We obtain analytical re- 
sults for bid-ask spread properties, Hurst plots and price increment correlation 
functions. 



1. Introduction 

While physicists involved in academic research on financial markets have focused 
on the analysis and modelling of time series such as stock market prices or 
indices |^|, |ll[ , limit order markets are addressed by economists in terms of market 
micro-structure pO[] . In such markets, traders can publish their wish to buy/sell a 
given quantity of stock at a given price, that is, place a limit order. Most of the time, 
there is no order of the opposite type at that price, and the limit order stays in the order 
book until it is cancelled or matched by another order of the opposite type at the same 
price. An order that causes an immediate transaction, as in the last example, is called 
market order. One of the most interesting aspects of these markets is the interplay 
between price dynamics and order dynamics. Whereas economists have studied and 
characterised many important aspects of these markets |2(J , such as why traders would 
like to use them after all ||, and produced models with stochastic order placement 
(see for instance || physicists prefer to consider orders as particles, thinking of 
size as mass and price as spatial position; the price evolution merely results from 
stochastic dynamics: it follows a kind of random walk, hence is not assumed to be 
fixed or bounded in contrast with other models p5| ^|, || . Data analysis performed in 
[jTo , 24 , |, [|, [35) shows that in the language of physics, limit order markets are out-of- 




equilibrium, that is, there are similar to a bucket of water heated by a campfire and 
filled by an irregular shower. There are two families of nonequilibrium particle models 
of such markets, one that assumes that orders diffuse on the price axis |l , [ 
and another one where orders can be placed, executed || and cancelled 1 1, |ll 
There is evidence from market data that orders do not diffuse on the price axis, but are 
placed, and then executed or cancelled with time dependent rates || . This begs for a 
non-equilibrium model made up of these three processes || |l2|, |J]. So far, all models 
have assumed constant rates. Here we show that this assumption is responsible for 
very unrealistic price behaviour, and that lifting this assumption leads to the correct 
behaviour. 
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It is well documented that stock and foreign exchange prices are over-diffusive 
over several months [^6|, [□]]: if the (log)-price p t followed a standard random walk, 
its increments (or returns) would be such that p T ~ Po ~ t h with H = 1/2; H is 
called the Hurst exponent. Generally speaking, H depends on r. For instance, real 
prices are reportedly over-diffusive (H > 1/2) for several months, and then tend to 
be diffusive (H = 1/2) for larger time lags [^6[ plj . Over-diffusion implies some sort 
of long-term memory. This behaviour is not captured in the current particle models. 
On the contrary, H = 1/4 for all times in models where no cancellation of orders is 
allowed 0, [li], |23|. It has been realised that order cancellation is needed in order to 
obtain a crossover to diffusive price for large time lags ||. On the other hand, over- 
diffusive behaviour is obtained when price trends are introduced, that is, when the 
price has a temporary tendency to rise or decrease, as in the model of Ref . , which 
has H = 2/3 for all times, when no order cancellation is allowed, and a crossover to 
H = 1/2 when it is. Here we generalise the notion of trends introduced in ]33| , and 
study their influence on the Hurst exponent and the return correlation function. 

Nonequilibrium particle models are notoriously hard to tackle analytically, and 
limit order market models are no exception. Therefore one ought to study first the 
simplest models of this kind, namely exclusion models p3|j , where at most one order 
is present at each tick. Recent work in nonequilibrium statistical mechanics shows 
that some of these particle models can be exactly solved [|3| (see ||l| for a review). 
Exclusion is naturally found when there are so few orders that it is unlikely to find 
two orders at the same tick, which is also the limiting case where Ref obtained 
analytical results concerning order densities in a class of models introduced in || [i~2| . 
Beyond this seductive mathematical aspect, exclusion models are arguably relevant 
for the modelling of limit order market. Indeed, priority of execution is given to orders 
first by price and then by chronological order of their deposition. This is particularly 
important near the best prices, when impatient traders are likely to be willing to take 
precedence over their colleagues. Therefore exclusion models of limit order markets 
capture at least some part of the dynamics near the bid-ask spread, which plays a 
major role in these markets, particularly when there is no market maker. Another 
interpretation of the exclusion condition is that it encodes the presence or absence of 
quotes at a given tick and neglects the actual number of shares. 

In short, the exclusion hypothesis is a way to reduce the dimensionality, hence 
the complexity, of limit order market models, while keeping their collective character. 
Here we transform the particle model found in Ref. Q by identifying processes in real 
markets into an exclusion model. 

2. Definition of the model 

Two types of orders, ask (^4) and bid (B) , also denoted by — and + and corresponding 
respectively to sell and buy orders, live on a semi-infinite line — the price whose unit 
is the tick. As limit orders are a way of storing demand and offer, it is clear that all bid 
orders are found on the left of the ask orders. At time t, there are ni(t) 6 {—1, 0, 1} 
orders on site i, where —1 denotes the presence of an ask order and 1 the presence of a 
bid order. We denote by b(t) the position of best bid price and a(t) the position of the 
best ask price. The last paid price p(t) is the position where a transaction last took 
place. All quantities of interest depend on these last three quantities. For instance, 
the bid-ask spread is the gap that separates the best bid from the best ask; here we 
denote it by s(i) = a(t) — b(t) — 1. We do not consider any mechanism which would 
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ensure the positivity of the price, since we aim to study normal price dynamics; we 
argue that special mechanisms arise in real markets when the price reaches very low 
prices.| We assume that all processes occur relative to the current best prices, that 
is, relative to aft) for asks and b(t) for bids. At each time step t, three different events 
can happen: a new order of each type (bid+ask) is placed, a market order of each 
type is placed, and some orders are cancelled. Mathematically, 

(i) with probability <5 a , a new ask order is placed into the bid-ask spread. Its price 
a(t) — A (A > 0) is randomly drawn according to some distribution probability 
P a (A). With probability 5' a , another ask order is deposited on the right of the 
current best price at a(t) + A', where A' > is randomly drawn from P^(A'). 
Idem for bids: an new order is placed in the spread with probability 8b at b(t) + A, 
where A is taken from P,(A). With probability S' b , a new order is placed at a 
lower price that b(t), that is, at b{t) — A' where A' is drawn from the probability 
distribution function (pdf) P b '(A'). There can be at most one order per site, 
hence, if there is already a particle at b(t) — A', or at a(t) + A', nothing happens. 

(ii) Each particle has the same fixed probability rj to be cancelled. 

(iii) With probability a a (t) a market order eats the best ask order located at a(t) and 
with probability ab{t) the bid order at b(t) is executed. 

(iv) Rates change according to predefined rules (see below). 

Deposition and market orders in particle limit order market models can be found 
in Ref. (23), cancellation was introduced in Ref. j(| and in Ref. in a different context; 
its main role is to ensure that the total number of orders does not diverge with time. 
Note that P(A) has a major influence on the bid-ask spread dynamics; the average 
typical deposition position conditional to the bid-ask spread depends linearly on the 
spread in real markets ||, which we denote by (A 2 \s(t)} 1 ^ 2 = cq + cis(£), where (.) 
denote a temporal average. This can be explained for instance by buy orders that are 
placed at b(t) + s(t)/2, or very close to the best ask order, e.g. at a(t) — I. Playing 
with co and Ci modifies the spread dynamics. For instance, if co > 0, ci = 0, there is 
no feed-back of the spread on the width of deposition; if c\ > 1, the contribution of 
cross-deposition to market orders cannot be neglected. Finally, if cq = and c\ > 0, 
the cross-deposition rate is constant. 

Several simplifying hypotheses will be made here: we will assume that all 
deposition pdfs are equal, that is, P a = p, = P' a = P fc ', and Gaussian. This contrasts 
with empirical facts that P'(A) is a power-law over two decades of ticks in the Paris 
Bourse (| m the London Stock and in the NASDAQ which implies that 

the probability of far deposition is non-ncgligiblc. In addition, the cancellation rate 
of a given order has a power-law dependence on its position relative to the best order, 
77(A) ~ A _l/ with v ~ 2 ^7j, which is consistent with the functional form of the 
lifetime distribution previously reported jl9| |): here, this dependence is neglected. 

In real markets, rates are not fixed, but vary wildly during any trading day. 
Moreover, we see no practical reasons why the rates should obey a balance equation 
such that the expected price drift is zero at all times, as it is assumed in all current 
models, except that of Ref. fl33"[ . On the contrary, heterogeneity of market participants 
is probably enough to argue that rates should be balanced only on average. As a first 

I There is technical reason to justify this assertion at least on the NASDAQ: a stock is unlisted from 
the NASDAQ if its price stays for too long under 1$ (see www.nasdaq.com). This is a motivation for 
the company to weigh on its stock price. 
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approximation, we assume here that they are independent from each other. As we shall 
see, this provides a universal mechanism for mid-term over-diffusive prices. We study 
several alternative simple rules that stipulate how to change the rates, depending on 
the quantity to be computed: 

(i) At each time step, with probability p, all rates are redrawn independently at 
random. The Hurst exponent is then very easy to compute. 

(ii) At each time step, each rate is redrawn with probability p. This is also appropriate 
for the computation of the Hurst exponent. 

(iii) Rates are linearly cross-correlated, that is, 



where o~(t) — ±1 is the only quantity that is changed with probability p at each 
time step. It defines the direction of trend. This setup makes it easy to compute 
the autocorrelation function and the properties of the bid-ask spread. 

(iv) Same as rule (iii), with the additional constraint that ao = a\ and So = 8\. In 
this case, market orders occur at one side, and order deposition at the other. 

(v) Same as rule (iv), but the direction of the trend a{t) is allowed to change only 
when the bid-ask spread is zero ]33| |. 

In rules (i-iv) the probability of rate change at each time step is constant. This 
sets the typical lifetime of a set of rates to 1/p, and that rates are exponentially cross- 
and auto-correlated, whereas this dependence should be algebraical |l6). Rule 
(v) relates trend changing to the actual properties of the bid and ask distribution. 
Other similar rules could also be studied, for instance by relating probabilistically the 
trend to volume imbalance between the two distributions. However, rule (v) is most 
appropriate for mathematical analysis. 

3. Bid-ask spread properties 

Here we derive basic results about the bid-ask spread. We assume that the width 
of order deposition does not depend on the bid-ask spread (ci = 0), and that the 
price of new orders new prices a(t) — 1 and b(t) + 1, that is, -P(A) = S<\,i where S 
is the Kronecker index. Finally, we neglect cancellation, whose role is less crucial in 
exclusion models, as the exclusion condition itself limits the total number of orders. 
The last two assumptions imply that after a transient fl36[ , the order density is uniform 
and equal to 1 for p < b(t) and p > a(t). 

As long as the rates are fixed, the best ask evolution is the following: P[a(t + 1) = 
a(t) + 1] = a a , and P[a(t + 1) = a(t) — 1] = S a , hence a(t) follows a biased random 
walk (RW) with drift a a — 5 a - The same applies to b(t), with a drift of 6b — cib 1 
hence the mid-price, defined as m(t) = (a(t) + b(t))/2, diffuses with an average drift 
(a a + 5f, — «b — 5 a )/2. The evolution of the spread s(t) is also a biased RW, with drift 
v = a a + at — 8 a — 5b < 0. Since s(t) > 0, the spread undergo a biased RW with 
reflective boundary at s = 0. Several properties of s can then be easily computed. For 
instance the pdf of the time during which s > is given by the first passage at s = 
of a particle starting at sq @1 



cv. a (t) = aa + Oi±a(t) ab(t) = ao — a\a{t) 
S a (t) =5 - 5xo-(t) S b (t) =5 + 6ia(t), 



(1) 

(2) 



(S + Vt) 

4Dt 



2 1 



F(t) 
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Figure 1. Trend duration pdf F(0,t), for rule (iv) (s = 1,5= 0.5, a = 0.20 
(squares) and a = 0.4& (circles)). 10 7 iterations. Continuous lines are theoretical 
predictions from Eq (p|). 



where D is the diffusion constant. In particular, when v — > 0_, F(0,t) — i~ 3 / 2 . 
For any v < 0, it is well known that the pdf P(s,t) of random walks with reflective 
boundary tends to a stationary state that is exponentially decreasing 



P(s) 



D 



exp(— \v\s/D) 



(4) 



and s = D/\v\, where ? denotes an average over P(s).y As shown by these equations, 
the spread properties only depend on v and D, provided that v and D remain constant 
long enough. 

Rules (iii), (iv) and (v) for rate change are such that the change of trend a(t) has 
no influence on the dynamics of s(t), as v and D are independent of o~(t), therefore 
the above equations are always valid in this case. In particular, rule (iv) makes the 
spread dynamics a one-step process for s > 0: 
dP(s) 



dt 



= [(1 - 6){1 - a) + ad]P(s) 
+ 6(1 - a)P(s + 1) + (1 - 



S)aP(s - 1) (5) 

Note that for s = 0, any deposition is in fact a cross-deposition, hence the effective 
annihilation rate a c g equals a + 8. Eq. (JsJ) shows again that the spread process is 
equivalent to a biased RW with reflective boundary at s = 0; its drift is v — a — 6 < 0, 
and its diffusion constant D = [a(l — a) + 8(1 — S)]/2. If order cancellation is allowed, 
the 0-th order correction to the above picture is 5 — > 6 — rj and a — > a + r/, hence 
v — > v — 2rj. 



4. Hurst exponent 

The Hurst exponent H depends in general of r in our model. Irrespective of the 
considered rule for the update of the rates, the price is over-diffusive (H > 1/2) for 

§ In continuous space and time, P(s, t) is indeed the solution of the diffusion equation Dd 2 P/ds 2 = 
-\v\dP/dt. 
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Figure 2. Hurst plot of the price increments for rule (i). . Market order rates 
drawn uniformly between and 0.1, deposition rate between and 1 (circles), 
and 0.5 (squares), p = 0.001, 10 6 iterations, no order cancellation. The dashed 
line represents normal random walk (H = 1/2) and dotted lines are obtained from 
Eq. <§). 

intermediate r and diffusive (H = 1/2) for large r. For rules (i)-(iv), the crossover 
between H > 1/2 and H — 1/2 occurs for r ~ 1/p. When rule (v) is applied, it occurs 
for t of the order J dtF(0, t)t. The small r behaviour depends on the rule and on the 
parameters, as discussed below. 

4.1. Rule (1) 

The simplest setup for computing the Hurst exponent H(t), is to change the rates 
according to rule (i). For a given set of rates a a , atb, S a , 5b, the average speed of the 
mid-price m(t) is v = a a + 8b — ctb — S a , Neglecting the diffusion of the mid-price, 
and retaining only its ballistic motion, one readily computes (r(t + r)r(t)) ~ (1 — p) T . 
Therefore, after a summation, 

{(p t+T - Pt f) = {C£ r t+v f) ~ r + \pr - 1 + (1 - pf] , (6) 

T^x p 

Figure || shows that Eq (jfy describes adequately the price behaviour for large r. The 
discrepancy observed for small r is probably due to bid-ask spread bouncing and 
other diffusive behaviours that we neglected. Note that when the deposition rates are 
close on average to the market order rates, a sub-diffusive behaviour appears for small 
times. This general behaviour is due to spread bouncing, as in this case, the spread 
is typically large. It appears for all rules where frequent spread bouncing is allowed, 
owing to the possibility of having simultaneous market orders at both sides, that is, 
rules (i), (ii) and (iii). 

4-2. Rule (ii) 

Rule (ii) is arguably more realistic than rule (i) . While it leads to essentially the same 
behaviour, the calculus is more complicated and of no special interest for our purpose. 
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Figure 3. Hurst plot for rule (iii), showing characteristic under-diffusive prices 
at short times. (8 = 0.5, a = 0.22, ai = 0.02, i) = 0, p = 0.001, 0.01 and 0.1, 
10 6 iterations) 



4.3. Rule (Hi) 

Fig. [3] shows that, similarly to rule (i) and (ii), if p is small enough, the price is 
under-diffusive, then over-diffusive, and tends to be diffusive at large times. Since the 
duration of the over-diffusive part is related to 1/p, when p is too large, the price 
cannot develop its over-diffusive behaviour. Note that the crossover from H > 1/2 to 
H = 1/2 has the same characteristic form as in Figs || and ^. 



4-4- Rule (iv) 

In the case of Rule (iv), the price follows a persistent RW p2\ , that is, the probability 
that r(t + 1) = r(t) is equal to I — p, where p is the probability of trend reversal. 
Therefore, the over-diffusive behaviour of our model can be qualitatively explained 
by mapping the price increments (±1) to an unidimensional (nearest neighbour) Ising 
spin model where the space plays the role of the time. A spin is a variable that can 
take two values —I and +1 (think of a price increment). The unidimensional Ising 
model consists in placing a given number of these on a line, and to assuming that spin 
at site i is only influenced by its neighbours, that is, the spins at sites i + 1 and i — 1 
(see |m| for more information). The mapping is possible, because the market model 
we consider is Markovian. Generalising known results for such models to our case 
where the price jumps over the spread when the trend reverses, whose effect turns out 
to be negligible for large r, we find when p is small and (s) C 1/p that for large t, 

, r(l-p) 



p(t)Y 



a 



P 



• tanh[rp/ (1 — p)] 



» 1 (7) 



which indicates a continuous crossover between H(t) — 1 for small r and H = 1/2 
for large t. However, the crossover appears visually to occur in a small time window. 
The mapping between the market and Ising models gives good results, even for small 
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Figure 4. Hurst plot for rule (iv) for a = 0.9, <5 = 1 (stars), a = 0.44, 
S = 0.5 (circles), a = 0.22, S = 0.25 (squares), and a = 0.08, <5 = 0.1 (triangles) 
(p = 0.001, r) = and 10 6 iterations in all datasets). 

4-5. Rule (v) 

This rule also gives rise to a persistent RW, hence to the same results as those of rule 
(iv). 

5. Price autocorrelation 

In this section, we study how the price increments autocorrelation function C(r) = 
(r(t + r)r(t)) depends on the rule of rate changes and on the chosen parameters. 

5.1. Rule (iv) 

Suppose that the probability of trend reversal is equal to p at each time step. 
Neglecting the effect of cross-deposition when s — 0, we find 

C(l) ~ C (l) = a 2 [(l - P f -p(l -p)(l - 6) 

-p*((s*) + (s)(l-6))}. (8) 

As soon as the influence of spread bouncing is greater than that of annihilation, 
C(l) < 0. This is to be related to one known possible cause for negative C(l), spread 
bouncing [^9[ [TTJ , although here it only happens whenever the trend a(t) changes. 
When P S (Q) is not negligible, Eq (||) must be corrected: a transaction takes place if a 
market order is placed (a), or a market order is deposited onto the best opposite order 
(S) (cross-deposition) . If both happen during a time-step, we assume that annihilation 
occurring before and after deposition are equally likely. We finally find 

C(1)~(1-P S (0))C (1) 

+ P,(0){(1 - p) 2 a[a + 6(1 - a/2)(l + 5/2)] 

+ p{l - p)[a 2 6/2 -a{a + S + a6 2 /A)} - p 2 a5 2 } (9) 
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Figure 5. P rice increment autocorrelation function C(l) 



subsection 5.1 



of the model of 
i, 5 = 0.4. r) = 0, 



(left panel: 8 = 0.5, a = 0.4; right panel: 5 = 
average over 50 runs, 10 s iterations; continuous lines are theoretical predictions 
of Eq (EH), and dashed lines are those of Eq. (H)). 



Figure ^| shows that numerical simulations are in good agreement with Eq (g) . Note 
that for this figure we used values of P s (0), (s) and (s 2 ) directly taken from the 
simulations, as the continuous approximation for P s breaks down for typically small 
spreads. In the right panel, the typical spreads is very small and Co(l) is not a good 
approximation to C(l). Eq (g) relates the price increments autocorrelation function 
to the average spread and the second moment of the spread, hence, can be seen as a 
generalisation of the seminal result of Ref [^9| , and is related other works J2l|, [l7| . 

5.2. Rule (Hi) 

Simultaneous transactions on both sides of the market are an additional source for 
short term negative return autocorrelation. Results of the previous subsection can 
be generalised to this case: the spread dynamics is described by a two-step process, 
which is again equivalent to a biased random walk with reflective boundary at s = 0, 
and C(l) can be calculated, but leads to expressions too long to be reported here. 

5.3. Rule (v) 

Ref. p3[ studies a particle exclusion market model where the order diffuse toward 
best prices and the rates are changed according to rule (v). According to the previous 
section, the trend duration pdf in this model is given by Eq (g) (see also Fig. 0) ; m 
addition, the rates are tuned so as to ensure that v = 0. Therefore the trend duration 
has no definite average, which explains intuitively why without order cancellation the 
price is always over-diffusive in this model; on the other hand, when order cancellation 
is allowed, the effective rates have to be corrected; this yields an effective |v| = 2r\ > 0, 
which is the origin of the crossover to diffusive prices observed in this case. Because 
this rule does not allow spread bouncing, the absolute value of price increments is 
or 1, and the sign of the increment stays constant for a duration F(l) that is given by 
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Eq (||). According to Ref. [|l5], the increment autocorrelation C(r) of such processes 
is the inverse Laplace transform of 



C(s) = I 



1 n 1-^) 



!(l)(l + F( S )) 



(10) 



where F(l) is the pdf of trend duration I. Unfortunately, C(s) cannot be inverted 
for F given by Eq (0). However, one can argue from this equation that for r small 
enough, C(r) ~ t -1 ^, following Ref. Bj. Another possibility is to consider a standard 
biased RW starting from sq and compute the probability that s(t) < 0, neglecting the 
effects of the reflective boundary, which leads to 



so 



11) 



(r(t)r(t + r)) = 1 - 2P[r(t) = -r(t + r)] ~ erf ^ _^ 

where the last approximation is valid as long as vt + s$ 3> 0, that is, as long as the 
spread is not likely to have reached the reflective boundary. When r is such that 
< 1, that is r < 2D/|u|, (r(t)r(t + r)) cx r" 1 / 2 , with a cutoff at i - 2D/\v\. 

Note that this result is exact if v = 0, and gives (r(t)r(t + r)) ~ r -1 / 2 for r ^> Sq/2D, 
in agreement with Ref |33|. The problem of this approach is that it gives positive 
autocorrelation for all times, whereas it is negative at short times in real markets p|. 

6. Price increment distribution 

All rules, except rule (v), lead to the same kind of price increment pdf, which is made 
up of one part that is due to consecutive market orders on the same side, and another 
that is due to spread bouncing (the price pdf produced by rule (v) lacks the spread 
bouncing part). This implies that when the rates are constant, the tails of P(r) are 
given by that of the spread pdf, which decreases exponentially in that case. In other 
models, found e.g. in Refs |^3|, [|, the pdf of the spread is fat-tailed, hence that of 
the prices as well, which is arguably artificial, since the mid-price increments m(t) are 
also fat tailed in real markets pl]| . 

It is possible to obtain a more refined pdf of the trends. Up to here we essentially 
considered constant rates, which lead to exponential tails. Time varying rates is the 
key to obtain fat tailed return pdf in these models: the mid-spread point m moves 
indeed at average speed ( s B- a B)+(aA-s A ) _ ^ + ai + § Q + 5i)<j. The rates are in 
principle not bounded, as it is easy to generalise our model to include multiple order 
depositions/annihilations during a time step. In this case, they indicate the average 
number of events during At = 1 of a Poisson process. Their distribution is related: 
about 10% of deposited orders are matched in the London Stock Exchange |L0| as 
well as in the ISLAND ECN (www.island.com), a subpart of the NASDAQ §], hence 
one can consider the deposition rate as proportional to the annihilation rate in a first 
approximation. P{r) has fat tails as soon as P(a) ~ a -7 . According to Ref. pi , 
this is the case, with an exponent of about 3.5. In addition, the autocorrelation of all 
the rates was found to be algebraically decaying which also shows up in that of 
the number of transaction in a given time interval [ lq] . We shall explore further this 
generalisation in a forthcoming publication. 
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Figure 6. Integrated average bid density (B)(A) versus A, the position relative 
to the best bid for r) = 0.01 and spread feed-back (S = 0.5 a = 0, co = 2, ei = 3). 
The dashed line has an exponent of -1. 



7. Average densities 

Two recent papers [0, ^] have focused on the average densities of orders relative to 
the best prices. In both cases, these were found to have power law tails. The model of 
Ref. H and its exclusion version, i.e. the model proposed here with spread-feedback, 
are able to reproduce this behaviour (see fig |J);|jj] (B(A)) is very well fitted with the 
functional shape (1 + ax) b with b ~ —2.03, very close to —2. This gives a market 
impact exponent of 1/2, which was measured in real markets, as reported Ref p8|. 



Note that the model of Ref. |12 also gives the same exponent: this model can be 
seen as a special case c± = l/v2, cq = of Ref. m, with a flat spatial probability 
distribution of deposition; in particular, spread-feedback is present (ci > 0), which 
seems to be necessary for obtaining densities with power-law tails. 

Interestingly, the functional shape of the average densities is the same as that of 
the price in these models. In fact, it is also the same as that of the spread, a general 
feature of market models with constant rates, as explained above. This last feature is 
probably very unrealistic, as the pdf of the spread is most likely exponential in real 
markets. Nevertheless, this can be easily explained. Suppose that the spread pdf has 
power tails, then order deposition pdf has fat tails, hence, according to the argument 
developed in ref Q|, the average densities have power-law tails. In addition, as before, 
when the rates are constant, the tails of the price pdf are due to spread bouncing, 
hence are also fat. In conclusion, in this type of models, the spread, the price and the 
densities have all the same type of pdf. But as argued above, fat tailed return should 
be obtained by modulating the rates. Note that the model of Ref. Q| produces power- 
law tails for the densities by considering power-law distributed relative positions for 
bulk deposition, which is consistent with empirical findings. 

|| One can even simplify further the exclusion model by replacing one of the two order distributions 
by a wall, still obtaining this result. 
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8. Conclusions 

The simplest particles models of limit order markets are exclusion models and provide 
an ideal framework for mathematical analysis. We argued that particle models of 
limit order market should include time-varying rates, and that the balance of rate 
only holds on average, and is a crucial feature of these markets, explaining mid-term 
over-diffusive price behaviour. We introduced several rules for time varying order 
deposition and market order rates and characterised the bid-ask spread properties, 
computed Hurst plots and correlation functions. The effects of order cancellation and 
more general probability distributions of rates will be addressed in the future. 

This work was supported by EPSRC under Oxford Condensed Matter Theory 
grant GR/M04426. We thank G. M. Schiitz for seminal discussions. 
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